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Force-Free  Magnetic  Fields, 
Curl  Eigenfunctions,  and  the 
Sphere  in  Transform  Space, 
with  Applications  to 
Fluid  Dynamics  and 
Electromagnetic  Theory 


1.  INTRODUCTION 

More  than  two  decades  ago  Moses1  developed  representations  for  eigenfunctions  of  the  curl  operator, 
applying  them  to  a  number  of  physical  problems  of  interest.  We  employ  these  functions  here  to 
construct  a  new  basis  for  the  representation  of  force-free  magnetic  fields  that  will  facilitate  their 
application  to  the  solution  of  geophysical  and  astrophysical  problems. 

The  concept  of  the  force-free  magnetic  field  (FFMF)  has  proven  very  useful  in  modeling  a  number  of 
physical  systems  of  interest  since  its  introduction  by  Lundquist2.  the  initial  development  by  Lust  and 
Schliiter3,  and  the  appearance  of  the  useful  mathematical  description  by  Chandrasekhar  and  Kendall4. 
While  physical  applications  have  ranged  from  stellar  astrophysics  to  controlled  thermonuclear  fusion, 
the  dominant  topic  of  papers  published  has  been  the  modeling  of  solar  magnetic  fields  and  their  effects 
on  the  upper  atmosphere  of  the  sun5.  The  magnetic  clouds  ejected  from  the  sun  w  hich  have  produced 
the  major  perturbations  to  the  Earth's  radiation  belts  during  the  satellite  era  seem  to  possess  FFMFs 
which  have  budded  from  the  solar  magnetic  field  (Burlaga  and  Lepping6,  Klein  and  Burlagu7. 
Lepping,  Jones,  and  Burlaga8).  The  present  work,  which  is  motivated  by  the  desire  to  describe  these 
solar  magnetic  clouds,  their  origins,  and  their  interactions  with  the  Earth's  magnetosphere,  introduces 
a  new  viewpoint  for  the  investigation  of  these  fields  and  indicates  its  broader  application  in  fluid 
dynamics  and  electromagnetics. 
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In  spite  of  the  numerous  papers  which  have  dealt  competently  with  various  parts  of  the  theory,  the 
mathematical  description  of  FFMFs  remains  incomplete  with  many  ramifications  unexplored.  (As  a 
cogent  example,  note  the  recent  appearance  of  a  modern  analysis  of  properties  of  the  curl  operator. 
Yoshida  and  Giga9.)  The  purpose  of  the  present  paper  is  threefold:  to  introduce  a  new  representation 
ot  these  fields  which  offers  both  practical  computational  advantages  and  a  change  in  theoretical 
view  point,  to  suggest  a  scheme  for  classifying  FFMFs.  and  to  indicate  the  broad  applicability  of  the 
new  representation  to  fluid  dynamics  and  electromagnetics. 


2.  THE  GENERAL  FORCE-FREE  MAGNETIC  FIELD 

The  concept  of  a  force-free  magnetic  field  arose  naturally  in  describing  the  equilibrium  states  of  a  non- 
resistive  plasma.  Considering  the  particular  limiting  case  in  which  the  hydrostatic  pressure  and  the 


Lorentz  force  are  both  zero 

V/i  =  J  x  B  =  0  .  ( 1 ) 

Lundquist2  realized  that  the  requirement  that  the  current  density  satisfy  the  static  Maxwell  equation 

V  x  B  =  n  J  (2) 

amounted  to  a  constraint  on  the  magnetic  field 

(V  x  B)  x  B  =0  (3) 

which  may  be  equivalently  stated 

(V  x  B)  =  a  B  .  (4) 

Thus,  for  the  magnetic  field  to  be  force-free,  the  magnetic  field  vorticity  must  be  parallel  to  the  field 


itself  with  a  proportionality  that  is  in  general  a  function  of  the  spatial  coordinates: 

a  =  a(x)  (5) 

The  type  of  functional  dependence  of  a  partitions  the  class  of  force-free  magnetic  fields  into  two  sub¬ 
classes:  variable  a  characterizes  the  so-called  nonlinear  subclass,  while  constant  a  characterizes  the 
other  subclass.  We  restrict  our  attention  here  to  this  latter  case,  as  have  most  previous  investigators. 

The  investigation  presented  here  is  mathematical:  the  results  obtained  are  thus  independent  of  the 
degree  to  which  the  FFMF  model  is  applicable  to  any  given  system  and  thus  applicable  to  all 
appropriate  systems.  See  Roberts10,  Chapters  1  and  4,  for  a  helpful  discussion  of  the  physical 
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idealizations  and  approximations  commonly  used. 


3.  THE  CONSTANT-crFORCE-FREE  MAGNETIC  FIELD 

3.1  The  Mathematical  Boundary  Value  Problem 

The  mathematical  problem  addressed  here  is  the  determination  of  the  magnetic  field.  B.  which 
satisfies  the  first  order  partial  differential  '  quation  system 

V  x  B  =  a  B 

(6) 

V  B  =  0 

(with  a  constant)  together  with  appropriate  boundary  conditions. 

A  thorough  discussion  of  the  boundary  value  problem  (BVP)  is  complex  and  lengthy,  and  would 
divert  us  from  the  main  point  of  the  report;  it  is  deferred  to  a  later  report.  The  question  of  the  proper 
posing  of  the  BVP  may  be  perceived  as  an  undercurrent  running  through  many  of  the  published 
papers,  especially  those  dealing  with  the  numerical  modeling  of  the  fields.  The  early  intuitive 
comments  of  Grad  and  Rubin1 1  have  frequently  been  used  to  justify  the  approaches  taken  by  those 
later  authors  who  have  confronted  the  problem  directly.  The  literature  here  is  extensive;  consult  the 
review  by  Sakurai12  for  orientation,  and  the  following  papers  for  a  small  sampling  of  the  variety  of 
approaches;  Saks13,  Bostrom14,  Semel15,  Aly16,  and  Yan  et  al.17.  The  series  of  papers  by 
Kress18-19’20  present  a  careful  formulation  of  certain  FFMF  BVPs  in  terms  of  integral  equations. 

The  constant  parameter  a  plays  a  central  role  in  the  development  below.  Although  it  can  be  viewed 
merely  as  a  scale  factor  in  physical  space  (as  is  evident  on  inspection  of  the  system  (6)),  the  choice  of 
this  scale  factor  has  significant  implications,  as  will  be  seen  below.  It  was  recognized  early  that  a  is 
intimately  involved  with  the  boundary  conditions  imposed  on  the  field  ( Lust  and  Schliiter3. 
Chandrasekhar21),  but  even  today  this  relationship  is  not  well  understood.  Were  we  to  address  this 
here  it  would  involve  us  deeply  with  the  boundary  value  problem.  In  keeping  with  the  decision  to  defer 
the  BVP  discussion,  we  do  not  follow  this  path  here. 

The  mathematical  basis  for  the  approach  taken  here  may  be  found  in  the  works  of  Ehrenpreis22  and 
Palamodov23,  in  which  Fourier  transform  methods  are  applied  to  the  solution  of  systems  of  partial 
differential  equations  with  constant  coefficients.  We  assume  that  the  general  solutions  derived  below 
can  be  fitted  to  appropriate  boundary  values  using  the  inverse  Fourier  transform,  such  as  exemplified 
by  Sneddon24.  Since  our  primary  concern  is  the  elucidation  of  a  new  approach  to  FFMFs.  no  direct 
application  of  this  technique  is  presented  here. 
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Remark  There  are  interesting  ties  between  the  current  approach  and  a  number  of  areas  that  are  worth 
mentioning  en passant.  The  BVP  as  treated  by  Chandrasekhar  and  Kendall4  invokes  the  second  order 
vector  Helmholtz  equation  obtained  by  applying  the  curl  to  the  first  of  Eqs.  (6)  above.  The  field  is  then 
obtained  by  constructing  a  vector  solution  from  solutions  to  the  scalar  Helmholtz  equation  in  analogy 
with  the  vector  waves  found  by  Hansen25  and  developed  by  Elsasser26.  The  sphere  in  transform  space 
is  as  fundamental  a  concept  for  the  scalar  Helmholtz  equation  as  it  is  for  the  curl  eigenfunction 
representation  (Miller27,  Newton2*).  Algebraic  group  operators  restricted  to  the  sphere  (the  transform 
images  of  differential  group  operators  in  x-space)  can  probably  be  used  to  investigate  the  separability 
of  the  FFMF  equation  in  various  coordinate  systems  in  a  manner  completely  analogous  to  that 
expounded  by  Miller27.  In  this  connection  the  FFMF  group  structure  derived  by  Yemets'  and 
Kovbasenko20  should  be  helpful.  While  a  number  of  investigators  have  implicitly  recognized  the  role 
of  the  sphere  in  transform  space  (e.g.,  Meyer30,  Dombre  et  al.31),  and  others  have  approached  the  use 
of  curl  eigenfunctions  (McLaughlin  and  Pironneau32,  Constantin  and  Majda33),  it  is  the  explicit  use  of 
the  Moses  curl  eigenfunctions,  with  their  orthogonality,  completeness,  and  associated  vector  Fourier 
transforms,  that  has  made  possible  the  analysis  described  here. 


3.2.  Curl  Eigenfunction  Representation  of  Force-Free  Magnetic  Fields 

The  system  (6)  reveals  that  an  FFMF  must  possess  vorticity  and  no  divergence:  it  is  a  solenoidal. 
rotational  field.  Following  Moses’  extension  of  the  Helmholtz  theorem,  such  a  field  may  be 
represented  as  a  sum  over  the  rotational  curl  eigenfunctions  (see  the  Appendix  for  a  summary  of 
properties  of  the  curl  eigenfunctions): 

B(x)  =  Jd  *k  %*(x  I  k)  /?A(k)  (7) 


where  the  prime  on  the  summation  index  indicates  the  omission  of  the  A  =  0  (irrotational) 
eigenfunction,  restricting  the  summation  to  the  values  A  =  ±  1.  Substitution  of  this  representation  into 
the  defining  equation  ( 1 )  and  making  use  of  the  eigenproperties  of  Xa.(x  '  k)  leads  immediately  to  the 
equation 


d  'k  (A  k 


a)  1  k)  Mk)  =  0  . 


(K) 


Taking  the  spatial  scalar  product  of  this  with  x^(x  I  k')  and  using  the  orthogonality  of  the  curl 
eigenfunctions  transforms  the  relation  to 
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d^k  (A  k  -  a)  5(k  -k’)  5U^  bx( k)  =  0 


(Ji  k'  -  a)  bit( k')  =  0 


(9) 


with  solution  (remembering  /i  =  l//j  for  =  ±1 ) 


vo 


-fy *»<*> 


(10) 


wliere  we  have  made  use  of  the  radial  delta  function  (see  Barton34),  k’  is  a  vector  drawn  to  a  point  on 
the  unit  sphere,  and  ,v;i(K')  is  a  function  yet  to  be  determined.  This  is  the  central  result  which  provides 
the  basis  for  the  rest  of  the  report.  There  are  two  immediate  consequences  of  this  relation: 

(1)  an  arbitrary  constant-a  force-free  magnetic  field  is  defined  entirely  in  terms  of  the  values 
of  its  transform  on  the  sphere  k'  =  alfi  =  a  / j : 

(2)  since  k’  must  be  positive  (as  the  magnitude  of  the  real  vector  k’),  only  one  eigenfunction 
contributes  to  the  field:  the  one  whose  sign  matches  the  sign  of  or.  n  =  sgn(a). 


Reconsidering  the  defining  equation  (4),  we  see  that  in  the  constant-a  case,  u  is  merely  a  (signed) 
reciprocal  length  unit,  and  the  subsequent  formulation  can  be  simplified  by  introducing  the  non- 
dimensional  variable 

^  =  I  a\ x  (11) 


which,  reverting  to  A  as  the  designator  of  the  curl  eigenvector  compouc"..,  changes  Eqs.  (6)  to  the 
form 


V'xB  =  A  B 

(12) 

V'  B  =0 

vvhere  the  prime  indicates  differentiation  with  respect  to  with  the  result  that  the  unit  sphere  in 
transform  space  assumes  a  dominant  position  in  the  theory.  The  expression  for  the  FFMF  now 
simplifies  to: 
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Bi(5) 


j*  k  2dk  dQ  ^  ^  U ^  sa(k)  X\(£>  ‘  K) 
Xi(£>  *)  'A(k) 

Ctt)  3/i  f  d£2  eIK  *>  QaOc)  *a(k). 


(1?) 


where  d  Q  =  sin  rj  dr]  dy/  denotes  the  surface  element  on  the  unit  sphere  in  terms  of  the  polar  angle 
t]  and  the  azimuthal  angle  y/  : 


k  x 


=  A  a  K  E 
I  a| 


(14) 


Remark.  It  has  been  known  for  some  time  that  the  general  solution  to  the  scalar  Helmholtz  equation 
may  be  east  in  terms  of  a  Fourier  integral  over  the  unit  sphere  in  transform  space;  see  Miller27,  as  well 
as  the  recent  detailed  investigation  of  Nev/ton28.  The  essential  difference  between  the  FFMF  and 
Helmholtz  solutions  is  the  appearance  of  the  Qa(k)  matrices  in  the  FFMF  transform  expressions. 


Thus  the  general  three-dimensional  constant-a  FFMF  is  defined  entirely  by  the  variation  of  the  scalar 
function  .vA(ic)  =  ,vA(r|,  yt)  in  conjunction  with  the  curl  eigenfunction  Xa(x  I  k  }  over  the  unit  sphere 
in  transform  space,  a  reduction  of  the  3-D  problem  in  physical  space  to  a  2-D  pro  demfor  a  function 
of  compact  support  in  transform  space 

An  additional  simplification  appears  on  examining  the  requirement  that  the  FFMFs  be  real.  For  this  to 
be  so.  the  field  must  be  equal  to  its  conjugate: 


Ba($>  =  B,*^) 


(15) 


or 


J  dQ  e  1  K  5  Qa<k)  ,va(k) 


J  dQ  e  '  K  *  Qa*(k)  V(k)  . 


(16) 


Inverting  k  on  the  right-hand  side,  multiplying  by  e  ‘  K  5  and  integrating  over  £  and  ic*  leads  to 


Qa(k')  sk(K')  =  Qx*(  k')  va*(  k')  , 


(17) 
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which,  making  use  of  the  inversion  relation  between  the  Qx(see  Appendix), 


Qa*(-k)  =  - 


KT,  -  /  A  JC2 " 
v  x'i  +  i  A  k2  , 


Qa(k) 


(18) 


leads  to  the  restriction  on 


*a(k)  = 


f  X'  l  -  /  A  <2  ' 
,KT,  +  I  A  K2j 


SX*(-K) 


(19) 


or.  equivalently. 


*a(-k)  = 


,  +  i  X  k2) 


(20) 


This  means  that  the  complex  FFMF  transform  need  only  be  defined  on  a  hemisphere  of  the  unit 
sphere,  since  the  transform  on  the  "lower"  half-space  is  defined  when  the  conjugate  function  is 
known  on  the  "upper"  half-space.  This  is  precisely  what  Moses  found  in  the  general  vector  field  case1. 

The  general  expression  for  the  FFMFs  may  now  be  simplified  by  using  the  above  relation  to  restrict 
(he  transform  to  a  hemisphere  in  transform  space.  Consider  the  general  expression: 

Ba($)  =  (2/r)  3/2  f  dGe,'**QA(K)5;l(K)  (21) 


For  ease  of  calculation,  we  orient  the  z-axis  along  defining  the  polar  angle  i)  as  the  angle  between 
k  and  £: 


cos  r|  s  k ■%!  £, 


(22) 


and  we  split  the  integration  over  the  sphere  into  that  over  the  +z  and  -z  hemispheres: 
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Ba(4)  =  (2 


sin  t]  dr\  e  '^C1,s  n  Qa(k)  sa(k) 


+  1  sin  Tj  d  7]  e  ‘  ^  cost,Qa(k)  sa(k)] 

* 


(2?) 


Now  invert  Kin  the  -r  integral,  setting 

k  =  —  k1  ,  r]  =  n  -  T]1.  \f/  =  n  +  y/  (  24) 

to  obtain 


P*  _  wfl 


Ba(4)  =  (2 


sin  rj  dr\  e  '^c<>sTI  Qa(k)  ja(k) 


»[} 

•  0  W  0 

+  sin  T]‘  d  i]'  e  ‘  ^  cos  f)'  Qa(-k’)  sa(~k’)  ] 


(25) 


and  make  use  of  the  reality  conditions  to  find 


Qa-Ki  aa(-k')  = 

^  K  j  + 

(  K  '  I  -  1  A  K  S  ^ 

Qa*(k')  . !  2 

,  K'l  -  i  A  K-';  , 

Vk'i  +  i  A  k'2  y 

,rA*(K')  =  Q  A*(k’)  sa*(k') 


(26) 


leading  to  the  expression 


=  (2/t) 


'"4 

*  0  *0 


sin  t]  di)  c'S tt,s  n  Qa(k)  ja(k) 


+  J*  sin  rj'  dr]'  e  '  £cos  H' Qa*(k')  ,va*(k')  J 


(27) 
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which  is  obviously  real,  since  the  second  integral  is  clearly  the  conjugate  of  the  first,  and  taken  only 
over  the  limits  of  the  upper  z -hemisphere. 

Thus  we  have  shown  that  a  g.  nu«l  constant- a  force-free  magnetic field  is  defined  completely  by  its 
variation  on  a  unit  hemisphere  in  Fourier  tranform  space. 


4.  APPLICATIONS  OF  THE  CURL  EIGENFUNCTION  REPRESENTATION 

We  now  demonstrate  that  the  representation  introduced  above  readily  simplifies  and  clarifies  certain 
FFMF  concepts  and  is  capable  of  generalized  application  to  fluid  dynamics  and  electromagnetics. 


4.1.  Vector  Potentials  for  the  Force-Free  Magnetic  Field 


The  current  density  and  the  magnetic  field  must  both  be  nonzero  w  ithin  the  region  considered  in  order 
for  an  FFMF  to  exist.  For  nonzero  a,  an  FFMF  must  be  described  by  a  vector  potential,  since 
otherwise  the  system  (6)  becomes  inconsistent.  In  the  curl  eigenfunction  representation,  the  vector 
potential  is  in  general  a  sum  of  two  rotational  modes  analogous  to  Eq.  (7)  (equivalent  to  adopting  the 
Lorentz  gauge  in  the  absence  of  an  electric  potential): 


A(x)  = 


d  k  Xa(x  I  k)  nA(k) 


(28) 


If  this  is  now  equated  to  an  FFMF  mode 


Ba<x)  =  V  x  A(x)  = 


?*F' 


k  k  Xa(x  I  k)  aA(k) 


(29) 


then  the  spatial  scalar  product  leads  to  the  expected  relation  between  the  field  and  vector  potential 
transfonns  : 


byik)  =  Xkax(  k)  =  aaA(k). 


(30) 


The  modes  must  match.  On  transforming  to  x-space,  we  find  the  corresponding  relation 


J 


BA(x)  =  Id  k  />A(k)  JCa<x  I  k)  - 


a  J"  d^k  a 


*(k)XA(x  I  k)  =  qAa(x) 


(31) 


As  Moses  found  in  the  general  case,  there  is  no  mixing  of  modes;  i.e.,  a  magnetic  field  described  by  an 
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eigenfunction  with  A  =  +1  must  be  derived  from  a  vector  potential  also  described  by  an  eigenfunction 
with  A  =  +1.  For  an  FFMF  this  becomes  a  more  stringent  result:  the  only  curl  eigenfunction 
component  needed  to  describe  the  vector  potential  of  an  FFMF  is  the  one  whose  eigenvalue  matches 
that  of  the  magnetic  field.  After  all,  a  constant  s  FFMF  is  fundamentally  a  curl  eigenfunction,  and 
vice  versa. 


4.2.  Magnetic  Helicity  of  the  Force-Free  Magnetic  Field 

Several  papers  by  Chandrasekhar  and  Woltjer35  and  by  Woltjer36-3738  (see  also  the  recent  work  of 
Laurence  and  Avellaneda31')  were  devoted  to  exploring  the  energy  and  stability  of  an  FFMF  by 
examining  the  behavior  of  a  pseudoscalar  quantity  which  we  denote  by  M: 


1 


M  —  I  d  x  A(x)  B(x) 


n2\ 


where  the  integral  is  taken  over  the  space  occupied  by  the  field.  This  quantity,  termed  the  magnetic 
helicity  by  Moffatt40,  has  a  fundamental  importance  for  any  magnetic  field  (see  the  discussion  by 
Ranada41  and  references  therein,  and  Berger  and  Field42),  but  takes  on  a  simple  form  for  a  constant- 
aFFMF  in  the  curl  eigenfunction  representation.where  it  takes  the  form 


Mx  =  J  Aa(x)  Ba(x)  =  a  J 


d^x  Aa(x)'? 


(33) 


or.  inserting  the  transformed  expression,  we  find: 


Mx  =  a  J*  d*k  |  a-A( k )  j' 


(?4) 


Thus  the  magnetic  helicity  for  the  constant-a  FFMF  is  proportional  to  the  mean  square  magnitude  of 
the  vector  potential  taken  over  r-space  or  to  the  mean  square  vector  potential  transform  taken  over  k- 
space  with  a  proportionality  constant  equal  to  a  (known  in  differential  geometry  as  the  abnormality 
of  the  vector  field  B  (or  A .  since  they  satisfy  the  same  equation  for  the  given  gauge  choice):  see 
Ericksen4'  and  Truesdell44).  The  abnormality,  a, may  be  defined  to  be  the  normalized  magneto 
helicity,  M A: 


ot  -  M  x  - 


J 

S 


d  x  A  Ax)  B a(x  ) 


d  x  A  A(x)  A  a(x) 


(35) 
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The  magnetic  helicity  is  but  one  element  of  an  infinite  helicity  hierarchy  whose  individual  terms  can  be 
very  simply  expressed  in  terms  of  the  curl  eigenfunctions  in  the  constant-a  case.  We  consider  this 
hierarchy,  couched  in  terms  of  fluid  dynamics,  in  the  following  section. 

4.3.  Trkalian  Motion  in  Fluid  Dynamics 

It  is  well  known  that  the  general  force-free  magnetic  field  is  equivalent  to  the  general  Beltrami  How  in 
fluid  dynamics,  while  the  constant-a  FFMF  is  equivalent  to  the  Trkalian  subset  of  the  Beltrami  flow 
(Beltrami45,  Trkal4(\  Truesdell44,  Bjprgum47.  Bjqrgum  and  Godal48,  Godal49).  Although  the 
development  above  is  couched  in  terms  of  the  magnetic  field  problem,  everything  discussed  there  can 
be  interpreted  in  fluid  dynamical  terms  with  the  velocity  of  the  fluid  acting  as  the  counterpart  of  the 
magnetic  field.  In  particular,  the  Trkalian  problem  can  be  formulated  entirely  in  terms  of  the  behavior 
of  the  velocity  transform  on  the  unit  sphere  in  transform  space.  Note,  however,  that  the  stability 
analysis  for  the  fluid  case  must  be  distinguished  from  the  analogous  magnetohydrodynamic  case,  as 
pointed  out  by  Moffatt  -50-51. 

In  this  section  we  consider  a  sequence  of  integrals  that  are  second-degree  functions  of  the  fluid  velocity 
and  which  we  call,  for  want  of  a  better  term,  the  helicity  hierarchy.  In  keeping  with  the  restriction  of 
this  paper  to  constant-a  fields,  the  hierarchy  is  considered  for  Trkalian  fluids,  though  the  formal 
application  to  general  Beltrami  flows  is  straightforward  and  will  be  addressed  in  a  later  report. 

In  analogy  to  the  definition  of  helicity,  we  introduce  the  helicity  hierarchy  as  the  following  set  of 
integrals  of  the  scalar  product  of  the  various  interated  curls  of  the  fluid  velocity: 

Hmn  =  f  d  \r  |  Vx  ]"'u(x)  •  [  Vx  jnu(x)  (36) 


where  ( m ,  n)  are  taken  to  be  (positive)  integers.  The  (kinematic)  helicity  is  immediately  recognizable 
as  H0 j.  Three  other  low  order  elements  of  this  sequence  occur  more  or  less  frequently  in  the 
literature:  the  (incompressible)  kinetic  energy  per  unit  mass.  H0 0;  the  enstrophy.  // s  j .  and  the 
superhelicity,  //12  (Hide52).  A  given  member  of  the  set,  Hmn ,  can  be  described  in  statistical  fluid 
dynamic  terms  as  the  correlation  of  the  mth  iterated  curl  of  the  velocity  with  the  nth  interated  curl, 
and  we  note  that  the  matrix  of  the  Hmn  is  symmetric.  We  now  consider  the  form  which  this  hierarchy 
takes  for  Trkalian  fluids. 

A  Trkalian  fluid  possesses  a  velocity,  u,  that  satisfies  the  system  (6)  above: 
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V  x  u  =  a  u 

V  u  =  0 


(37) 


hence  the  analysis  of  Section  3  goes  through  unchanged  and  we  can  write 


«a(x)  = 


J 


dQ  Sx{Q)  Xa(£  1  K) 


(38) 


In  view  of  the  properties  of  the  curl  eigenfunctions,  we  can  immediately  write  down  the  expression  for 
the  nth  interated  curl  of  the  velocity: 


Vx  \muk(x) 


h  km  Ux(k)Xxi*  'k) 
dQ  Sx(Q)  X&  I  K) 


(39) 


enabling  us  to  evaluate  the  helicity  hierarchy  defined  above  as  the  sequence  of  integrals 


H 


inn 


dQ'  Sa(K’) 


xl^1  K)  Xa(^!  k') 


(40) 


This  is  a  remarkable  result.  It  says  that  for  a  Trkalian  velocity  field  the  elements  of  the  set  of  fluid 
correlation  functions  defined  by  Eq.  (36)  are  all  reducible  to  H qq,  the  kinetic  energy  per  unit  mass  and 
a  power  of  the  velocity  field's  abnormality  (or  normalized  helicity).  The  magnitude  of  Hmn  grows  or 
diminishes  with  the  order  ( rn+n )  according  as  or  is  greater  or  less  than  unity.  The  extreme  simplicity 
of  the  helicity  hierarchy  for  this  case  emphasizes  both  the  very  special  nature  of  this  type  of  field  and 
its  intimate  connection  to  the  concept  of  helicity. 

Normalizing  the  HIWI  analogously  to  the  magnetic  helicity  of  the  previous  section,  we  find 
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(41) 


lj  _  nm+n 
nm„  -  CL 

so  that  the  first  four  named  normalized  correlations  become: 
kinetic  energy:  Hoq  =  1; 
kinematic  helicity:  Hq\  =  a; 
enstropliy:  H \\  =  a2; 
superhelicity:  Hn  =  a3. 

Thus  the  energy  and  enstrophy  are  positive  definite,  while  the  kinematic  helicity  and  superhelicitj 
adopt  the  sign  of  a.  For  negative  a,  the  elements  of  odd  index  ( m  +  n )  are  all  negative  and  the 
even  elements  are  all  positive,  while  for  positive  a.  all  elements  are  positive. 

The  analogous  magnetic  helicity  hierarchy  is  obtained  by  replacing  the  velocity  by  the  magnetic  field, 
though  the  indices  differ  for  the  unnormalized  elements  because  the  analogue  of  the  fluid  velocity  in  the 
definition  of  magnetic  helicity  is  the  vector  potential,  not  the  magnetic  field.  Of  course,  a  cross  helicity 
hierarchy  can  be  defined  in  terms  of  the  convolutions  of  the  fluid  and  magnetic  transforms,  but  we 
shall  not  pursue  that  here. 


4.4  Bivector  Form  of  the  Maxwell  Equations  in  the  FFMF  Representation 

The  representation  developed  above  for  the  force-free  magnetic  field  can  be  applied  directly  to 
electromagnetic  theory  to  produce  a  new  and  potentially  very  useful  viewpoint,  which  is  introduced 
here  and  will  be  developed  in  a  subsequent  report. 

At  the  turn  of  the  century  Siiberstein53  introduced  a  bivector  formulation  of  the  Maxwell  equations  by 
defining  a  complex  field  vector  with  real  and  imaginary  parts  given  by  the  electric  field,  E,  and 
magnetic  field,  B,  resp.  The  Maxwell  equation  system  then  takes  a  form  that  is  easily  treated  by  the 
FFMF  formalism.  Following  Stratton54,  but  changing  notation  slightly  to  avoid  a  conflict  w  ith  that 
used  above,  we  introduce  the  field  bivector 

G  s  B  +  /  E/c  (42) 

where  we  use  rationalized  MKS  units  and  c  is  understood  to  refer  to  the  speed  of  light  in  the  medium: 
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The  Maxwell  equations  then  take  the  form 


VxG  +  l  — 
c  at 


p. j ,  v  g  =  imp. 


(44) 


where  (p,  J)  are  the  charge  and  current  source  densities.  In  the  absence  of  sources  they  become 


VxG  + i =  0 ,  V  G  =  0 , 
c  at 


(45) 


while  for  harmonic  time  dependence  (read  equivalently:  Fourier  analysis  in  time)  their  formal  similarity 
with  the  FFMF  equations  is  striking: 

VxG  =  aemG  ,  V  (»  =  0  .  (46) 

where  we  have  defined  the  electromagnetic  helicity  as  a,„  s  The  analysis  of  Section  3  is  fully 
applicable  here,  and  we  find  that  the  three-dimensional  electromagnetic  wave  problem  is  completely 
defined  on  the  sphere  in  transform  space  for  which 

k  =  A  aem  (47) 

w  ith  the  fundamental  difference  from  the  FFMF  problem  that  the  full  sphere  must  be  utilized  because 
we  cannot  require  the  field  to  be  real;  the  real  and  imaginary  parts  of  the  vector  G  are  independent, 
being  the  magnetic  and  electric  field  vectors.  Turning  this  result  around,  we  see  that  the  real  and 
imaginary  parts  of  an  FFMF  may  he  interpreted  as  the  magnetic  and  electric  fields  of  a  corresponding 
electromagnetic  problem. 

This  result  can  be  generalized  immediately  (by  an  appropriate  redefinition  of  aem  and  G)  to  encompass 
chiral  media  characterized  by  the  Drude-Bom-Fedorov  constitutive  relations,  simplifying  significantly 
the  equation  system  obtained  recently  by  Lakhtakia55  (see  also  Lakhtakia.  Varadan  and  Varadan56). 

Moreover,  the  full  BVP  for  the  electromagnetic  field  should  be  accessible  from  this  treatment  through 
appropriate  Green's  function  techniques.  If  further  work  confirms  the  utility  of  this  approach,  it  could 
have  a  significant  impact  on  both  the  theory  and  practice  of  electromagnetics. 

We  leave  the  mathematical  problem  at  this  point,  having  established  the  formal  connection  between  the 
more  genera!  electromagnetic  theory  and  the  force-free  magnetic  field  formulation,  the  latter  clearly  a 
subset  of  the  former.  Although  the  bivector  description  of  the  electromagnetic  field  has  been  known  for 
some  time  (Weber57,  Silberstein55,  Bateman58,  Stratton54),  the  FFMF  connection  has  not  been 


previously  pointed  out. 


5.  CLASSIFICATION  OF  CONST  ANT-a  FORCE-FREE  MAGNETIC  FIELDS 

Having  found  the  general  approach  to  constant-a  force-free  magnetic  fields  outlined  above,  it  is 
natural  to  ask  whether  it  can  be  used  to  elucidate  general  properties  of  these  fields.  In  particular,  does  it 
suggest  a  basis  for  classifying  the  various  types  of  FFMFs?  Such  a  classification  can  be  approached  in 
various  ways,  depending  in  part  on  the  motivation  of  the  classifier.  For  example,  one  type  of 
classification  might  be  based  on  the  generalized  multipole  expansion  of  Lust  and  Schliiter3  (see  also 
Lust,  Schluter  &  Trefftz59),  arranging  fields  into  classes  according  to  the  parametric  dependence  of 
their  multipole  coefficients.  By  analogy  with  electrostatic  multipoles,  such  an  approach  would  be 
expected  to  be  useful  in  describing  FFMFs  that  are  concentrated  in  limited  regions  of  space.  In  contrast 
to  such  a  scheme,  which  relies  on  the  characteristics  of  the  fields  in  coordinate  space,  the  classification 
we  sketch  here  is  based  on  the  dimensionality,  symmetry,  and  complexity  of  the  FFMF  transforms  on 
the  sphere  in  transform  space.  Whether  this  provides  a  useful  classification  of  types  remains  to  be 
proven  by  subsequent  development,  but  the  insight  gained  in  the  attempt  already  seems  quite  valuable. 

Viewed  at  the  most  general  level,  any  function  or  distribution,  defined  on  a  set  of  points  on  the 
transform  sphere,  which  possesses  a  curl  transform  (as  we  denote  the  Fourier  transform  containing  the 
complex  vectors  QA)  satisfying  the  defining  system  (6)  is  a  valid  FFMF.  We  partition  the  class  of 
FFMF  transforms  first  of  all  according  to  the  dimension  of  the  supporting  set  into  point,  curve  and 
area  transforms,  or  sets  of  dimension  0,  1,2,  as  well  as  possibly  fractal  dimension.  The  further 
partitioning  of  the  transforms  within  these  groups  is  done  in  terms  of  complexity,  aided  by 
considerations  of  group  symmetry.  Only  the  barest  outline  of  this  program  can  given  here:  details  will 
be  given  elsewhere.  To  provide  a  hint  at  the  results,  examples  of  the  simplest  fields  in  each  category 
are  exhibited  with  brief  indications  of  their  properties. 


5.1  Point  Transform  Fields 

The  simplest  point  transform,  corresponding  to  a  single  point,  is  a  delta  function  characterized  by  a 
given  amplitude  and  unit  vector  drawn  to  a  point  on  the  transform  sphere: 

n<fl>  =  %  -  50  ^  -  n")  **  *•'  (48) 

sin  r\  sin  i] 

leading  to  a  field  expression 
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(49) 


Ba(5)  =  (2 /r)  3/2  I  Qa(k)  s0*K  ~  K()) 

J  sin  ?) 

=  (2jt)  V2s0e  '  K#^  Qa(Ko) 

whose  real  or  imaginary  pan  can  be  recognized  as  the  standing  wave  that  is  the  well-known  simplest 
FFMF,  discussed  in  Moffatt40  (with  the  wave  vector  icq  oriented  along  the  positive  z-axis).  An 
ordered  hierarchy  of  point  fields  may  be  constructed  using  the  fundamental  expression  (49)  by  simply 
summing  the  fields  corresponding  to  different  point  constellations.  In  so  doing,  some  surprises  are 
found.  While  the  lines  of  force  of  the  single  point  field  (49)  are  straight  lines,  the  lines  of  force  of  a  2- 
point  quadrature  field  are  already  rather  complex,  and  those  of  a  3-point  quadrature  field  are  in  part 
chaotic  (and  can  be  related  to  the  ABC  field  suggested  as  a  prototypical  turbulent  velocity  field  in  fluid 
dynamics;  see  Amol’d60  and  Dombre31).  This  is  consistent  with  the  viewpoint  introduced  by  Sweet61, 
and  developed  by  others  (e.g.,  Cary  and  Littlejohn62),  that  the  magnetic  field  lines  of  force  are 
equivalent  to  the  phase-space  trajectories  of  a  dynamical  system  in  which  the  increase  in  complexity  of 
the  dynamical  system  leads  to  a  progressive  loss  of  integrals  of  the  motion  and  the  appearance  of  chaos 
(  Zaslavsky  et  al6i).  What  new  features  in  complexity  may  be  introduced  in  higher  order  members  of 
this  hierarchy  remain  to  be  discovered. 


5.2  Curve  Transform  Fields 

For  this  group  of  fields  the  transform  is  defined  on  a  continuous  curve  on  the  sphere  so  that  our 
reliance  on  distributions  can  give  way  to  the  use  of  ordinary  functions.  At  the  same  time  we  must 
actively  take  into  consideration  the  variation  of  the  transform  amplitude  on  the  defining  curves.  It 
seems  natural  to  consider  first  those  transforms  that  have  constant  amplitude  on  the  curve,  then  to 
expand  the  treatment  to  encompass  amplitude  variation.  As  an  example,  consider  the  simplest  family  of 
curve  transforms,  those  for  which  the  amplitude  is  constant  on  a  circle  of  constant  polar  angle  i ]0: 


sjv(fl) 


S(  T)-t]0) 

•y0  . . . 

sin  rj 


(50) 


After  evaluating  the  integrals,  the  complex  field  takes  the  form 

B*(x)  =  'A0  e,z//-'  (  (cos  t)(1+  l)Jn(p/Lp)ax-  (cos  tj0  -  l)J2(p/LP)('  ,2A*aA* 

2/rl/2 

-V2  X  I  sin  ?!„./,  (p/Lp)elXv k  ] 


(51) 


where  the  constant  vector  aA  is  defined  as 
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the  length  scales  are 


Lp 

L: 


_ 1 _ 

a  sin  7]0 


_ 1_  ... 

a  cos  7)0 


(53) 


and  (p,  <p,  z)  are  circular  cylindrical  coordinates  in  x-space.  and  the  Jn  are  the  usual  Bessel 
functions.  In  the  limit  of  rj0  — »  0.  Lp  — >  <»,  Lz  — »  1/a,  we  recover  the  point  field  found  above 
which  is  periodic  in  z  but  has  no  p  or  <p  dependence,  while  as  r)0  — >  nil,  Lp  — >  1/a,  L. 

°°,  and  we  find  a  completely  new  field  without  z  dependence  but  which  is  periodic  in  <p  and  quasi- 
periodic  in  p. 


5.3  Area  Transform  Fields 


These  fields  correspond  to  transforms  defined  on  2-dimensional  subsets  of  the  sphere.  The  simplest  is 
constant  on  the  sphere  and  is  equivalent  to  the  Fourier  transform  of  the  Q^: 

Ba(x)  =  (2n)  1/2  bw  [jo(ar)  +  ij\(ar)  Fi(cos  0)]  a*  -  ji  1/2  A  /  P]\(cos  0)  e  "t  <pj\(ar)  k 

(54) 

oo 

-  \A  X  H)"  +  l)jn(ar)  —  P* (cos  9)  e  ,2k(pax* 

V  n  =  2  («  +  2)! 

Here  r  represents  the  spherical  radius  and  the  jn  are  the  spherical  Bessel  functions,  reflecting  the 
essential  spherical  nature  of  this  field  (and  the  essential  spherical  nature  of  the  curl  eigenvector 
component  QA). 

To  handle  other  area  fields,  it  seems  expeditious  to  expand  the  transform  in  spherical  harmonics,  since 
any  function  defined  on  a  sphere  can  be  so  represented  (Courant  and  Hilbert64).  Other  functions  are 
also  available,  such  as  vector  spherical  harmonics  (Moses1)  and  Heun  functions  (Kalnins.  Miller.  & 
Tratnik65)  which  could  prove  more  useful  to  represent  functions  with  specific  symmetry  properties. 


5.4  Fractal  Transform  Fields 


The  previous  three  sections  have  exhibited  FFMFs  whose  defining  sets  on  the  transform  sphere  have 
integral  dimensions  of  0.  1  and  2.  Other  sets  possessing  non-integral  dimensions  are  easily  imagined 
within  the  context  of  modem  fractal  geometry.  Mandelbrot's  definition  of  fractal  is  "a  set  for  which  the 
Hausdorff-Besicovitch  dimension  strictly  exceeds  the  topological  dimension1'  (Mandelbrot66;  see  also 
Ed  gab’’7).  Fractal  support  sets  have  recently  been  introduced  for  dynamo  magnetic  field  models  (Finn 
and  Ott6*).  Whether  such  sets  have  an  important  role  to  play  in  the  further  development  of  FFMF 
theory  and  practice  remains  to  be  seen.  Their  construction  and  curl  transform  seem  non-trivial  and  are 
not  attempted  here,  but  the  possibility  of  their  existence  must  be  explicitly  acknowledged. 


6.  SUMMARY,  SPECULATIONS,  AND  FUTURE  DIRECTIONS 

6. 1  Summary 


By  introducing  the  Moses  curl  eigenfunctions  to  describe  constant-a  force-free  magnetic  fields,  we 
have  been  able  to  show  that  such  fields  are  defined  entirely  by  the  value  of  their  curl  transform  on  the 
unit  hemisphere  in  transform  space.  This  change  of  viewpoint  enables  an  orderly  exploration  and 
classification  of  their  properties,  an  introduction  to  which  was  briefly  sketched,  exhibiting  the  simplest 
fields  defined  on  0.  1,  and  2-dimensional  sets  on  the  transform  sphere.  The  new  representation  was 
applied  to  the  determination  of  several  quantities  of  interest:  (1)  the  magnetic  vector  potential,  showing 
no  mode  mixing;  (2)  the  magnetic  helicity,  demonstrating  that  the  parameter  a  can  be  interpreted  as  a 
normalized  magnetic  helicity;  (3)  the  helicity  hierarchy  (correlation  sequence)  for  a  Trkalian  fluid, 
proving  that  the  members  of  the  hierarchy  are  defined  entirely  in  terms  of  the  mean  kinetic  energy  and 
powers  of  the  parameter  a  Applying  the  representation  to  electromagnetics,  we  showed  that  the 
bivector  formulation  of  Maxwell’s  equations  reduces  to  a  force-free  magnetic  field  problem  in  which 
the  unit  sphere  in  transform  space  plays  an  essential  role,  opening  the  way  to  a  new'  approach  to 
electromagnetic  theory  and  applications. 


6.2  Speculations 

The  geometrical  viewpoint  employed  above  suggests  further  immediate  developments  and  applications 
of  the  curl  transform  representation,  several  of  which  are  being  actively  pursued  (see  §6.3).  We 
mention  briefly  here  connections  to  two  current  problems  that  are  more  speculative  in  nature. 

Turbulence;  Lagrangian,  Trkalian  and  magnetic.  The  form  of  the  helicity  hierarchy  derived 
above  undoubtedly  has  significant  implications  for  Trkalian  turbulence  (if  it  exists!),  making  the 
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determimtion  of  a  from  the  boundary  conditions  a  h  gh  priority  item.  While  a  connection  has  not 
been  attempted  here  between  the  helicity  hierarchy  and  direct  or  inverse  cascades  of  energy,  ensrrophy, 
heiicity  or  any  of  the  other  fluid  parameters,  it  bears  looking  inio  (cf.  Moffatt69  and  references  therein, 
especially:  Andre  and  Lesieur70;  Kraichnan71;  Frisch  et  al.72).  The  corresponding  effects  in 
magnetohydrodynamics  and  their  influence  on  the  generation,  maintenance  and  decay  of  dynamo  fields 
through  the  a-effect  of  Steenbeck,  Krause  and  Radler75  should  also  prove  interesting,  if  not  crucial 
(see  also:  Krause  and  Radler74).  The  time  dependence  of  Trkalian  flows,  if  similar  to  the 
corresponding  FFMFs,  may  be  simply  exponential  (Lundqutst2.  Chandrasekhar  and  Kendall4. 
Moffatt40;  see  also  Marris  and  Wang75,  ,.nd  Marris70-77).  If  more  complex  time  dependences  can  be 
found,  a  substantial  insight  into  the  cascades  and  other  characteristics  of  turbulence  and  dy  name  theory 
may  result.  The  question  of  stability  of  FFMFs  bears  on  this  and  has  generated  a  long  trail  of  papers  in 
the  literature.  A  substantial  fraction  of  these  deals  with  the  generation  of  current  sheets,  largely  in  the 
context  of  solar  flare  mechanisms  (sampling:  Syrov  at. skii78,  and  a  long  series  of  papers  by  Parker,  see 
Parker79  for  earlier  references).  Do  similar  localized  velocity  sheets  or  jets  exist  in  Trkalian/Beltrami 
flows  and  do  they  contribute  in  an  essential  way  to  turbulent  mechanisms? 

Laboratory  simulation  of  astrophysical  magnetic  fields.  The  most  difficult  part  of 
simulating  many  astrophysical  processes  in  the  laboratory  is  getting  the  scaling  of  the  essentia! 
parameters  correct.  If  this  is  not  done,  the  physical  processes  which  can  be  supported  and  studied  in 
the  laboratory  environment  may  bear  little  relation  to  those  of  the  real  universe.  Force-free  magnetic 
fields  may  turn  out  to  be  easier  to  model  in  the  laboratory  than  expected  because  of  the  dependence  on 
the  single  parameter  a,  which  is  itself  a  length  scale.  The  fixing  of  a  through  the  specification  of 
appropriate  boundary  conditions  will  probably  be  the  key  to  such  modeling.  W'e  note  -hat  transient 
force-free  magnetic  field  configurations  have  been  reported  recently  in  the  laboratory  by  Steaze!  and 
Urrutia80. 


6.3  Future  Directions 

The  current  work  provides  the  foundation  for  a  number  of  applications,  first  being  the  investigation 
and  classification  of  the  FFMFs  touched  on  in  Section  5  and  their  relation  to  problems  of  current 
interest,  especially  the  description  of  solar  magnetic  clouds  and  their  interaction  w  ith  the  Earth's 
magnetosphere.  A  complete  solution  to  the  FFMF  problem  can  only  be  given  when  the  influence  of 
boundary  conditions  on  the  field  form  is  knowm,  so  that  an  investigation  of  the  full  boundary  value 
problem  must  be  high  on  the  agenda.  The  enlargement  of  the  fluid  dynamic  investigation  to  establish 
connections  between  the  Trkalian  helicity  hierarchy  and  Lagrangian/Trkalian  turbulence  holds  out  hope 
of  near-term  progress.  Reconsidering  some  of  the  classic  electromagnetic  problems,  such  as  the  proper 
radiation  condition,  multipole  expansions,  and  the  description  of  polarization,  in  terms  of  the  FFMF 
formulation  of  the  bivector  Maxwell  equations  has  the  potential  of  providing  new  insights  to  old 
problems.  On  a  more  distant  time  scale,  the  application  of  the  present  formulation  to  the  dual  fluid 
velocity-magnetic  field  description  involved  in  the  dynamo  problem  would  seem  to  present  advantages 
due  to  the  natural  and  explicit  representation  of  the  essential  effects  of  vorticity  and  helicity. 
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Appendix 


Summary  of  Properties  of  the  Moses  Curl  Eigenfunctions 

Moses1  introduced  the  transform  space  vectors  defined  by 


Qo  (k) 
QaOO 


k 

k 

A_ 

vf 


k\  (k\  +  i  A  £2) 


-  1 


k{k  +  A.'j) 


A 2  (A  1  +  i  A  k  1 ) 
k(k  +  k$) 


-  i  A , 


A  t  +  i  A  A'  2 

K 


( A.l ) 


which  have  the  following  useful  properties: 


Qa*  (k)  •  Q/i  (k)  =  Sin  , 

X  Qa*  (k)  Qa>  (k)  =  Sij, 

A 

k  (k)  =  0  for  A  =  ±  1  , 
kxQi(k)  =  -»AQA(k)  . 


Using  these  vectors  he  defined  the  vector  functions  Xa 

-V? 

XA(x|k)  =  (2iz)  e‘k  *Q*(k)  (A.3) 


which  he  was  able  to  show  are  eigenfunctions  of  the  curl  operator: 


Vx&(x|k)  =  k  X  xa  (xlk) , 

V-XA(xlk)  =  0  for  A  =  ±l,  (A. 4) 

-3/2 

V-  Xo  (x|k)  =  -  i  k  (2n)  c‘k  x 

and  which  satisfy  the  following  orthogonality  and  completeness  relations: 
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(A. 5) 


j  dx  xx*  (x|k)  Xm  (xlk'l  =  5x»  <5(k  -  k')  , 

X  (  dk  Xk*  (*lk)  %x,  ( x ' I k)  =  S,j  8  (x  -  x’)  . 


The  latter  relations  enabled  him  to  show  that  any  vector  function  u  could  be  expanded  as 


u  (x) 


^k  Xx  (x!k)  Ux(k) 


(A. 6) 


It  is  the  latter  relation  that  is  central  to  the  current  report.  Moses  original  paper  contains  several  other 
significant  results,  including  an  extension  to  the  Helmholtz  vector  function  decomposition  theorem. 


I 


28 


